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, Abstract 

In this paper, we discuss the potential importance of measuring 
the CMB anisotropy bispectrum. We develop a formalism for comput- 
ing the bispectrum and for measuring it from microwave background 
maps. As an example, we compute the bispectrum resulting from 
""o ■ the 2nd order Rees-Sciama effect, and find that is undetectable with 

q current and upcoming missions. 

^ ; 1 Introduction 

> ■ 

^ . Observations of microwave background fluctuations are a powerful probe of 



the physical conditions in the early universe. Most analyses of the microwave 
background focus on measuring and interpreting the two-point function. If 
the microwave background fluctuations were a purely Gaussian random field, 
then the two-point function would completely characterize the microwave 
background. The non-linear growth of structure produces non-Gaussian fluc- 
tuations. If detectable, these non-Gaussian fluctuations are a new window 
into the evolution and growth of structure. 

The three-point function and its analog, the bispectrum, are potential 
tools for detecting these non-linear effects. At present, there is no clear sig- 
nature of such modes. There was no detection of three point correlations in 
the COBE data, nor any topological signatures of non-Gaussianityfrj. While 
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there is a reported detection of non-Gaussianity in the COBE bispectrum[0, 
||. Our ability to detect these effects will soon improve by many orders of 
magnitude with the launch of the MAP satellite (§ and the PLANCK || satel- 
lite. These upcoming launches motivate our detailed study of the bispectrum 
and the three-point function. 

In this article, we first discuss why the bispectrum is interesting (§ |2|), 
develop a general form for calculating the bispectrum (§ |3|) and then discuss 
the particular case of those perturbations resulting from the Rees-Sciama 
effect (§ £|). While the Rees-Sciama bispectrum is not detectable, we show 
in a companion paper that the non-Gaussian effects produced at low redshift 
through a coupling of gravitational lensing and either the Sunyaev-Zel'dovich 
effect or the late time ISW effect produces a detectable signal. 



2 Why is the Bispectrum Interesting? 

The three-point function and its transform, the bispectrum, are potential 
tools for detecting non-Gaussianity in the microwave sky. Just as the mul- 
tipole spectrum is a useful statistic for studying the two-point correlation, 
the bispectrum is useful in studying the three-point function. To put this in 
context, we define the ai m as the expansion of the temperature anisotropics 
into spherical harmonics: 



= J dhT(h)Y; m (h) (i) 



Thus, since the three point function is defined as B(l, m, n) = T(l)T(m)T(n), 
we define the bispectrum as: 

R mim 2 m 3 _ /r>\ 
hhh ~ u hm 1 U'l 2 m2 U hrn 3 \* ) 

The universe is thought to be rotationally invariant, and hence, a more 
useful quantity is the angle-averaged bispectrum: 

tj _ \ " / W ^2 ^3 \ nmim2m 3 fn\ 

^JiJaJa - I m m2 m I a hhh V 6 ) 

Even if there was no non-Gaussianity in the initial potential fluctuations, 
non-linear physics will produce mode-mode couplings. We argue here that 
the bispectrum is a particularly powerful tool for detecting these couplings. 
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In inflation, quantum fluctuations produce Gaussian random phase vari- 
ations in the gravitational potential: 



*(*) = / !^o(k)exp(zk-x) (4) 

with ($ (k)$g(k')) = P(A;)5 (3) (k — k'), the power spectrum of potential fluc- 
tuations at the current epoch. While processes during inflation may produce 
some non-Gaussianity, this is small in most inflationary models |J. On the 
other hand, when the fluctuations reenter the horizon, there are many non- 
linear effects that alter these initial conditions. Near decoupling, non-linear 
hydrodynamical effects couple modes 0. After decoupling, non-linear gen- 
eral relativistic corrections also alter the microwave background spectrum 
[§, H as do non-linear coupling between density fluctuations [|TU|, [□]]. How- 



ever, since potential fluctuations in the early universe are small, ($ 2 ) ~ 10 9 , 
these non-linear effects are usually viewed as undetectable. For example, 



Mollerach et al.[0] and Munshi et al. [ 11 ] estimated the reduced bispectrum 



generated by non-linear evolution of long wavelength modes and showed that 
they are undetectable in the COBE maps. In a spatially flat universe, the 
microwave background temperature at the surface of last scatter can be lin- 
early related to the potential fluctuations at that epoch, giving us what we 
shall henceforth call the linear temperature approximation: 

T L (n) = <j>{r r ) J —e^ Mk)g{k ), (5) 

where r is the comoving (conformal) lookback time, 0(r) is a normaliza- 
tion constant for potential perturbations which is equal to unity today, and 
g(k) is the linear radiation transfer function and describes the relationship 
between potential fluctuations and temperature fluctuations at recombina- 
tion. For small k, g(k) = 1/3. For larger k, we need to evolve the coupled 
matter-radiation fluid. Codes like CMBFAST compute g (k) for different 
cosmological models. 

Expanding the linear part of the temperature perturbation into spherical 
harmonics: 

«L = #Tr) J j^dn^Mmm^) . (6) 
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Since we will later wish to integrate over the sky, it will be useful to ex- 
pand out the exponential into spherical harmonics as well using the Rayleigh 
expansion: 

e ik ' r = 47r5:^(A;r)^(r)yi m (k). (7) 

Im 

Thus, the linear approximation of the CMB fluctuations are given by: 

at = 4vr0(T r )H)</ ^Mk)Y^ m m(kr r )g(k) (8) 

In addition to the terms proportional to $o(k) (linear), there are a number 
of terms in the temperature fluctuations which are proportional to $o(k) 2 
(nonlinear). These corrections arise from gravitational couplings, gravita- 
tional lensing, and radiation effects. These terms generically add a correction 
to the microwave background: 

T NL (h) = /°W 2 (n)/(nr,r) (9) 
= I rfT02(T) / (C)3 ( ^ k ) 3 $ o(k 1 )$o(k 2 )e-^^)- fi V(k 1 , k 2 , r) 

We may again take the spherical harmonic transform of the non-linear 
temperature perturbation to calculate the coefficient, . We will find it 
convenient to expand the coupling term, /(k 1 ,k 2 ,T) into Legendre polyno- 
mials: 

4-7T 

f(k u k 2 , r) = £ m, h, r)P,(k r k 2 ) = ]T ^—rfiih, h, r)Yi m (k 1 )F^(k 2 ) 

(10) 

Since the non-linear term will only couple to itself if k x = k 2 , we can safely 
ignore this effect, and instead concentrate on the case in which a bispectrum 
is produced by: 

TDhhh _ n L L NL* _i_ L L NL* , L L NL* (||\ 
L, mim 2 m3 ~~ Uj hmi Uj hm2 Uj hm i ' u 'l 2 m2 Uj hm s u 'hm 1 ' Uj hm i Uj hmi Uj hm2 \ LL ) 

We will only examine the first term in this derivation, and add in the other 
two permutation in the end. 
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Looking at a^*, we find: 



a ?£ = |^/^ 2 W/^i^2$o(ki)^(k 2 ) (12) 



(2vr) e 

X E jfl(h,k 2 ,T)ji,(k 1 T)j l »(k 2 T) 

WV'mm'm" 



x | dndkxdk.y^^ywtf^ 

Substituting this expression and equation (|8|) into equation (|TT|), we get: 



af imi af 2m2 aS = g^(r r ) / dr0 2 (r) / / J^fe^ (13) 

x $o(k)$ (k')$S(ki)$o(k 2 )^(A;)( ? (A:') 
x S 7^r—Tfi( k hk2,T)ji'(k 1 T)ji // (k 2 T)j h (kT r )ji 2 (k , T r ) 

IV l" mm' m" M + 1 

dndkiy / 3 m 3(n)y / / m /(n)y /w (n)y; m ,(ki)y / ^(k 1 )y;* mi (k) 

rfk 2 y m (k 2 )y;, m „(k 2 )y /2 * m2 (k) 

The integrals over k and k' give two possible permutations, <5(k — ki) and 
5(k — k 2 ). This is patently symmetric, and thus, we will chose the former 
and multiply the integral by 2. After integrating over k and k', we find that 
we have 3 sets of 3 spherical harmonics. This is precisely the form of the Q 
coefficients discussed in the appendix. Thus, integrating over n, k 1; and k 2 , 
summing over m,m',m" , and convolving with the Wigner 3-j symbol over 
mi,m 2 , rriz as in equation (^), we find that the first of three terms in the 
angle averaged bispectrum is: 

B hL = ^( T r) J ^(t) J dhdktfkZPikjPfaMkdgfa) (14) 

x £(2Z' + 1)(2*" + l)Q l ivrfi(ki, ^r^ihr^ihr^ihr^ihTr) 
ll'l" 



3 Measuring the Bispectrum 

Measuring this effect from observations is quite straightforward. If there is 
complete sky coverage and no sky cuts, the the angle averaged bispectrum 
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can be computed by rewriting the Wigner 3-j symbols as integrals over the 
sky, 



h h h 




B, 



El h h h 
I 

mi,m2,rri3 \ 



h h h 
mi n^2 



C^/imi ^?2"i2 



(4tt) ; 



\ (2Z 1 + 1)(2Z 2 + 1)(2Z 3 + 1) 

x E / rf q^im 1 (q)^ 2m2 (q)^ 3m3 (q) 

mi,m2,m3 

dqe h (q)e h (q)e l3 (q) (15) 



where the sky map is averaged over rings centered around point q: 



'2Z + 1 

47T 



(16) 



Since the bispectrum signal is rather weak, the dominant source of noise is 
the cosmic variance of the dominant Gaussian signal. For a Gaussian random 
field, {a hmi a hm2 ai 3m3 ) = 0. However, its variance, ((a limi a hm2 a hm3 ) 2 ) = 



c h c h c h f° r h 7^ ^2, h 7^ hi h 7^ ^3 [|T2|| . When we include detector noise, then 

{{a hmi a hm2 a hm . A f) = (c h + a^wf 2 ) (c h + crlw h 2 ) (c h + alw h 2 ) (17) 

where <To is the detector noise and wi is the experimental window function. 
Similarly, for a Gaussian field, 

( B hhh) = 

{B hhh B Vx i^ = (c h + alw h 2 ) (c h + (>lw h 2 ) (c h + a 2 w h 2 ) (18) 

x (4*2 + + <fi§ + <K + S ™ + &) 

for h ^ h, h 7^ h, h 7^ h- Because of this symmetry, we restrict ourselves to 
bispectrum terms where h < h < h- 
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4 The Bispectrum from the Rees-Sciama Ef- 
fect 



Several groups |I| [TT| have discussed the possibility of detecting the non- 
linear signature of the Rees-Sciama effect in the COBE data. They concluded 
that the signal was undetectable on COBE scales. In their analysis, they 
focus only on the diagonal terms in the bispectrum Here, we compute the 
predicted signal on smaller angular scales and compute all of the bispectrum 
terms. 

The second order Rees-Sciama effect arises from the non-linear growth of 
density fluctuations [pT0f| : 

d_ 



T NL (h) 



dT^ NL (hT,r) 



(19) 



2 / dr 



d 



d 3 k 



dr J (2tt) : 



$ JVi (k,r)e 



tk-nr 



and give the second-order potential perturbation as: 



(to 



84P 



(d 3 k') 

1>)3 



$o(k-k')$ (k / ) [3k 2 k' 2 + 7k 2 k ■ k' - 10(k ■ k') 

(20) 



l\2 



This second -order approximation is valid for k < lh Mpc" 1 . |[L3|| . 

Since we expect much of the signal to come from high redshift contribu- 
tions, and since in a flat, Q m = 1 universe, <p(r) = 0, we will take <f>(r) is 
constant for the remainder of this derivation. Additionally, we will redefine 
our variables of integration: ki = k — k' and k2 = k'. This yields: 

d 3 kx d 3 k 2 



T 



NL, 



n 



dr<p\r) 



Jo 



T 



21 k x 



(2tt) 3 (2tt) 3 
lOklkj + 7(kl 



$o(k 1 )$ (k 2 )e i(kl+k2)iiT 

k\)k x k 2 ^ + 4k\k\ii 2 



(21) 



where // = ki • k%. 

The second part of the equation is the coupling function, f(ki, k 2 , r). The 
vectors in the denominator make it difficult to decompose it into Legendre 
polynomials, however, this may be simplified by noting that: 



oc 



(22) 
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Thus, taking the Laplacian of a^, we find: 

Kl + 1)^ = / dhj dr^r) J ^^ $0 ( kl )$ (k 2 ) e ^+^y^(n) 

x /(k 1 ,k 2 ,r)(k 1 + k 2 )V (23) 

We now define a new coupling function: 

/(kx, k 2 , r) = /(kx, k 2 , r)(ki + k 2 ) 2 r 2 , (24) 

which can be easily substituted into the definition of the bispectrum, with 
B hhh => h(h + l)B h i 2 i 3 and /, =>- /,. 

For the Rees-Sciama effect, this is decomposed into only the I — 0,1 — 1, 
and / = 2 Legendre coefficients: 

. 2(r -r)r 2 



f (k 1 ,k 2 ,T) = 
fo(k 1 ,k 2 ,T) = 
f 2 (k 1 ,k 2 ,r) = 



3 

2(r - r)r 2 

3 ' 
2(to - r)r 2 





(25) 


\{k\+kl)) 


(26) 


Yi klk2 ) 


(27) 
(28) 



so, 



h(h+ l)B hhh = ^ 2 (r r ) j 'dr<f> 2 (r)(r -r)r 2 (29) 
x / dhdkikfklPikjPik^gik^gik^j^hT^ji^Tr) 



x Y,3i'{k l T)j vl {k 2 r){2l' + \){2l" + 1) 

VI" 



X 



2^Qoi'n' 3 kik 2 + -Qip/' 3 (^i + & 2 ) + -^Q 2 i»i?k\k 2 



As discussed in the appendix, £ n „(2Z' + l)(2Z"+l)Q^ 3 = (2Zi + l)(2Z 2 + 
l)<5o/il2 3 ' anc ^ thus, if / is small (in this case, it only takes on the values 0,1,2), 
and hi'i" is a function which varies slowly with I' and we may make the 
approximation: 

+ 1)(2Z" + l)Q|^ 3 ^p ~ (2k + 1)(2Z 2 + l)QSh l2h (30) 
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Under this assumption, the bispectrum reduces to: 

A(j) 2 (T r )I lll2h 



l 3 (k + l)B hl2l3 = 3^/2 V ^ 2(r)(r °- T)r2 (31) 
x J dk l dk2klklP(k 1 )P(k 2 )g{k l )g(k 2 ) 
x 3h (kir r )ji 2 (hr r )ji 2 (k 1 r)j h (k 2 r) (h + k 2 ) 2 



where 



4 w , = V(2/i + l)(2/ 2 + l)(2Z 3 + l)( o o o) (32) 



In order to deal with the spherical Bessel functions, we will first note that 
in the asymptotic limit, if x > I, ji(x) ~ sin(a; — nl/2)/x. The lower limit of 
integration will be taken to be: 

r ^= max (!'l) (33) 

Likewise, the lower limit of integration on k\ is l\jr T and k 2 is l 2 /r r . In this 
limit, we can expand out the four spherical Bessel functions as: 

n n 1/ , x7 -u c °s \(ki — A; 2 )r r )l cos \(ki — k 2 )r) 

MhTrfafaTrfaihT^faT) ~ -(-l)' 1+ ' a ^ ~ 

12 T 

(34) 

plus additional terms which oscillate much more quickly. 
We can now approximate the the time integral: 

/TV 
drr 2 (r - T)4> 2 (T)j h (k 1 T r )j h (k 2 T r )ji 2 (k 1 T)j h (k 2 T) (35) 
min 

^ ^-l^oosKfe-feK] r dr(ro-r)^(r) ro8[ ^"i 2)r] 

^ n Sjh-h) , )h+l 2( w v 

~~ 4 kfk 2 T 2 K'minJK'Q 'mm) 

where we have approximated 0(t)(to — t) as a slowly varying function which 
is maximally weighted near r = r min , and sin(x)/x is approximated as a 
delta function. 
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Thus, the bispectrum can be approximated as: 



f 37T 3 / 2 Z 3 (/ 3 + 1) Ji m / T 

(36) 

where we have defined l m = max(li, h)- 

In equation (H), we defined g(fc) in terms of o^. Defining q = (a^ m a{£), 
we find: 

Q = — / k 2 dkP(k)jf(kr r )g 2 (k) (37) 

For large values the spherical Bessel functions are very close to zero for 
kr r < I and since k 2 P(k) is a monotonically decreasing function for all values 
of k, we may approximate jf(x) = tt5(1 — x)/2. With this approximation: 

* ^ f »W (38) 

Substituting this into the previous expression yields: 

B hL - o L j TTTiS r dkk 2 P{k)c kTr {2kT r + l)<p\r mm ){T -T mm ) 
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Q0 (r r l m /l)(ro ~ TV WO 



3vr 3 / 2 r r 4 / 3 (/ 3 

r&U (39) 



1 ^iM 3 



37T3/2 r 4/ 3 (/3 + l) 

In Figure p] we have plotted the coefficient for flat cosmologies where 
Vt m = 0.1, 0.3, and 1.0. 

Will we be able to detect this effect? We can estimate this by determining 
whether in a Q m = 0.3, Q\ = 0.7 universe, we can use the MAP data to reject 
the hypothesis that there is no non-Gaussianity in the CMB maps: 

X 2 = E 1^4 (40) 

h,h,h \ h,h,hl 

In Figure |2| we show the amount of "information" gained by increasing the 
maximum index (/ 3 in our notation). Note that since all of the signal occurs 
around I < 600, and the MAP and PLANCK satellites are cosmic variance 
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dominated up to I ~ 600 and / ~ 1000 respectively, both experiments should 
measure the effect with approximately the same sensitivity. 

The integral of Figure |2| gives us the total x 2 between our fiducial model 
and some test model. In Figure |3|, we plot the total % 2 under the assumption 
that the "true" universe is a flat Q m = 0.3, and the test models remain flat 
and vary Q m . Even in the extreme case in which we assume Q m = 0.3 and 
the true universe has Q m = 1.0, the \ 2 for PLANCK is only ~ 0.02. 



5 Discussion 

Most previous discussions of the microwave background have focused on mul- 
tipole spectrum of the microwave background: 



1 



21 + 



- a im a *m (41) 



Because of azimuthal symmetry, the expectation value of (af m ) is independent 
of m. If the microwave background fluctuations were purely a Gaussian 
random field, then the multipole spectrum measures all of the statistical 
properties of the microwave background. 

In this paper, we have developed formalism for computing the bispectrum 
and for measuring it from microwave background maps. As an example, we 
computed the Rees-Sciama bispectrum and found that is undetectable with 
current and forthcoming missions. In a companion paper we will work 



out the bispectrum for nonlinear contributions from the coupling of gravita- 
tional lensing and low-redshift temperature perturbations. There, we show 
that MAP will be able to detect the SZ-lensing bispectrum and PLANCK 
will be able to detect the ISW-lensing bispectrum. 

Much work remains to be done on the bispectrum. We do not have a ro- 
bust statistical technique for computing the bispectrum from a million pixel 
map with sky cuts and spatially variable noise, nor do we understand how 
foregrounds will contaminate measurements of the bispectrum. For exam- 
ple, with MAP, radio sources are expected to produce marginally detectable 



skewness[15]. There are also other physical effects that may have observable 



signatures in the bispectrum. 
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A Wigner Symbol Identities 

In this appendix, we derive some relevant properties for the integral, 



ll'l" 




m u m 2 m 3 \ " b ± " l 2 " L 3 J 
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where cos a — I ■ m, cos (3 = rh ■ n and cos 7 = I ■ h and 

' I I" h \ ( V I h 



# 3 = J(47r)»(2/ 1 + l)(2/ 2 + 1)(2Z 3 + 1) 











I" V h 

(43) 



We can use the definition of the Wigner 6-j symbol [jig] to rewrite the integral 

as, 



r)hhh 
Wll'l" 



rhhh J ^1 k Z 3 

\ V I" I 

h k h 

V I" I 



-1) 



i+i'+i" 



E 



mi,m2r?i3 



h k k 

17li m 2 ^3 



■1) 



h+h+h 



rhhh 
1 ll'l" 



:-i) 



i+i'+i" 



(44) 



If I" = 0, then the only non-zero term is 



VM3O 



(2Z 2 + 1) (4tt)3 (l x l 2 l 3 
\ (2/ 1 + l)(2Z 3 + l) ^ 



(45) 



If I" = 1, Q\lri 3 is zero unless I — Zi ± 1 and Z' = Z 3 ± 1. Thus, there are only 
four terms in the sum in equation (27). 

We can evaluate these terms by noting that 



h-l l 2 Z3-I 




(/l+*2 + *3 + l)(*l-k + k) / h k h 
\ (h+h + h)(h-l2 + h-l) V 



h h-l 1 





(2Z 1 + 1)(2Z 1 -1) 



(46) 
(47) 



and 



Zi Z2 Z3 
Z3-I 1 Zx-1 



\ 



(k + h + h + l)(/i + / 2 + h)(h -k + h- l)(/i - + 



(2Zi + 1)2Z 1 (2Z 1 - 1)(2Z 3 + 1)2Z 3 (2Z 3 - 1) 

(48) 

Thus we can reduce the product of Wigner symbols to a simpler form: 



2(2Z 1 -l)(2Z 3 -l) ^ hhQ 



(49) 
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Similarly, 



I lJ2 i 3 (/l+/ 2 -Z 3 +l)Hl+/ 2 + /3) n M 2 /3 

ii+i,i3-i,i 2(2Zi + 3)(2Z 3 -l) V,1,s0 ' 



n iii 2 «3 _ (-^1 + ^2 + ^3 + + / 2 - h) nh l 2h 

- 2(2/i _ 1)(2/3 + g) Vl!i 3 > 

and 



hhh (h + l 2 + h + 2)(h-l 2 + l 3 + l) hl2h 

h+1,/3+1,1 2(2/ 1 + 3)(2/ 3 + 3) ^ lZ3 °' 



Note that 

£Q|^ 3 (2Z + 1)(2Z' + 1) = (2k + l)(2/ 3 + 
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Figure 1: The coefficients 6/ m for various cosmologies, as defined in the text. 
Note that for low Q models, the signal is significantly larger than Q m = 1 
models at all 1. 
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Figure 2: The amount of information gained by increasing / 3 , as given by 
X 2 = E((B hh i 3 {n m )) - (B hh i 3 (a m = 0.3))) 2 /(B? ihh ), where we have used 
the PLANCK detection sensitivity. Since there is very little contribution 
above 1% = 700, the MAP and PLANCK missions will be equally sensitive to 
this effect. Our fiducial model here and throughout is Q m = 0.3, Qa = 0.7. 
The shape of this spectrum comes from several competing effects. At large /, 
more modes are included in the sum, providing a larger signal. However, the 
signal per mode becomes increasingly weak at a logarithmic rate for high 1, as 
suggested by the plot of the coefficients in Figure 1. Finally, the features in 
the spectrum above are caused by variations in the expected noise. Troughs 
correspond to peaks in the CMB spectrum. 
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0.2 0.4 0.6 0.8 1 
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Figure 3: The total value of x 2 (fi m , Vt rn = 0.3) for variations in Q m for both 
PLANCK and MAP. All values of fl m produce signals which are consistent 
with Q m = 0.3. 
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